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■ The Equation of State (EOS) for asymmetric nuclear matter is discussed 
starting from a phenomenological hadronic field theory of Serot-Walecka type 

I including exchange terms. In a model with self interactions of the scalar 

' sigma-meson we show that the Fock terms naturally lead to isospin effects 

■ in the nuclear EOS. These effects are quite large and dominate over the 
contribution due to isovector mesons. We obtain a potential symmetry term of 



o 
o 



"stiff' type, i.e. increasing with baryon density and an interesting behaviour 
of neutron/proton effective masses of relevance for transport properties of 
asymmetric dense matter. 



^ ! Phenomenological hadronic field theories (Quantum Hadrodynamics, QHD) are widely used 
in dense nuclear matter studies since relativistic effects are expected to increase with baryon 
density 0]. In most of the previous works on the subject, the Relativistic Mean Field (RMF) 
^ ' approximation of QHD has been followed. In the RMF the meson fields are treated as 
classical fields and consequently a Hartree reduction of one body density matrices is used. 

Although the model has driven a large amount of results on relativistic effects in nuclear 
structure and dynamics [0-0, the lack of exchange terms has implied some non satisfying 
features of the theory and some efforts have been done to try to cure this problem |17|-[10|. 
In the RMF theory each meson field is introduced, with appropriated readjusted couplings, 
just to describe the dynamics of a corresponding degree of freedom, without mixing due to 
many-body effects. Neutral a and uj mesons are in charge of saturation properties, isospin 
effects are carried by isovector 6 [ao(980)] and p mesons and finally spin terms are coming 
from pseudoscalar vr and t] fields. In a sense the model represents a straightforward extension 
of the One-Boson-Exchange (OBE) description of nucleon-nucleon scattering. 

The aim of this letter is to introduce explicit many-body effects just evaluating exchange 
term contributions. We will get qualitative new features of equilibrium (EOS) and dynam- 
ical properties of asymmetric nuclear matter. In particular a new density dependence of 
the symmetry term is expected, at variance with the simple linear increase predicted by the 
p— exchange mechanism in the Hartree scheme. 
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As we already know from non-relativistic effective interactions, like the Skyrme forces, 
Fock terms play an essential role in symmetry breakings and consequent mixing of different 
degrees of freedom. Similar effects are expected here. In particular, in the context of the 
QHD model, essential properties of nuclear matter come mostly from the two neutral strong 
meson fields. Hence it is important to evaluate the Fock contribution associated with these 
fields. 

We will focus our attention on isospin contributions to the nuclear EOS, symmetry 
term and neutron/proton effective masses, and on the relativistic transport equation for 
asymmetric nuclear matter. 

We start from a QHD — II model where the nucleons are coupled to neutral scalar 
a and vector u mesons and to the isovector p meson. Self-interaction terms of the cr-field 



were originally introduced for renormalization reasons |]TT],|T2| and can also be considered 
as a way to parametrize the density dependence of the NN force. Actually they are also 
describing medium effects essential to reproduce important properties (compressibility and 
nucleon effective mass) of nuclear matter around saturation density. The Lagrangian density 
for this model is given by: 

C = ^[lf.{id^ - gvV^ - 9pB^ ■r)-{M- gs<P)]4j + ^{d,<j)d^<j) - m^V') 

~4>' - - Iw^uW^" + \mv^V,V^ - ^ V ■ + \rnp''B, ■ (1) 

where ly^^(x) = 9'^V^(x) - d^'V^^x) and L'^^(x) = (9'^i3^(x) - ^''B^'{x) . 

Here il){x) generally denotes the fermionic field, and V'^(x) represent neutral scalar 
and vector boson fields, respectively. B^{x) is the charged vector field and r denotes the 
isospin matrices. 

From previous equation one can derive the field equations and the canonical energy- 
momentum tensor [|I|. 

In our approach we will perform the many-body calculations in the quantum phase space 
introducing the Wigner transform of the one-body density matrix for the fermion field. This 
method has two main advantages, the use of physical quantities and the direct derivation of 
dynamical transport equations. The one-particle Wigner function is defined as: 

X - —) :> , 



[F{x,p)U = 7^ / d'Re-'^-'' <: + 



here a and (3 are indices for intrinsic degrees of freedom of the fermionic field (spin and 
isospin). The brackets denote statistical averaging and the double dots denote normal or- 
dering. The equation of motion for the Wigner function can be derived from the Dirac field 



equation by using standard procedures (see e.g. Refs. |]T3|,|T3), it reads: 

'-d^[YF{x,p)U+p^[YF{x,p)U - MF^p{x,p) 
(2^ / <■ M^MMx-Wx^) :> 

+937:^, f d'Re-'P-'' <: ^p{x+)Mx-)<P{x^) ■> 



(27r)4 . 

'^^ (2^ / <■ M^MMx-)r . B,{x^) :>= , (2) 
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with x+ = X + ^ and x^ = x — ^. 

In order to take into account the contribution of exchange terms in a manageable way 
we assume, as a basic approximation, that in the equations of motion for the meson fields 
the terms containing derivatives can be neglected with respect to the mass terms. Therefore 
the meson field operators are directly connected to the operators of the nucleon scalar and 
current densities: 

+ A$2 + 5$3 = iIj{x)iIj{x) , 

Vi'ix) = fv^{x)-ff'ip{x) , 

mx) = fpi,{x)^^T^{x) , (3) 

where fs = {gs/msf, fv = {gv/mv?^ fp = {gp/rripf and $ = gs(j), V^' = gyV, A = a/g^, 
B = b/g^, = gpB^. After substituting in Eq.(2) these expressions for the meson field 
operators, we obtain an equation which contains only nucleon field operators. 

The present approximation implies that retardation and finite size effects in the exchange 
of mesons between nucleons are neglected. However, we are concerned with a semiclassical 
description of nuclear dynamics, so that the nuclear medium is supposed to be in states 
for which the nucleon scalar and current densities are smooth functions of the space-time 
coordinates. Therefore, because of the small Compton wave-lengths of the mesons cr, uj and 
p, the assumptions expressed by Eq.s(3) are quite reasonable. For light mesons such as pions 
this approximation is not justified. The inclusion of self-interaction terms of the pionic field 
requires a different approximation scheme. However, it has been shown that the inclusion 
of pions does not change qualitatively the description of nuclear matter but, as expected, in 
very low density regions 0. 

An attempt to include exchange terms in the QHD approach was previously performed 
without self-interaction terms for the a field [^], with results not satisfying due to the 
inadequacy of the model. Within the same model, a relativistic kinetic equation with self- 
consistent mean field has been derived in Ref. |15| taking into account the exchange terms. 



Here we evaluate the effects in a more physical model with self-interacting higher order a 
terms. 

In order to evaluate the nuclear Equation of State the quantity of interest is the statistical 
average of the canonical energy-momentum density tensor. According to our approximation, 
where terms containing the derivatives of the meson fields are neglected, the tensor takes 
the form: 

T,uix) = '-^P{x)^,d,^{x) + [Ui$) - lv,V^/fv - ^B, ■ B'/fp]g,,, (4) 

where g^u is the diagonal metric tensor and U{^) = / fs + A/?) $^ + 5/4 <l>^. Following 
the treatment of the Fock terms in non-linear QHD introduced in Ref. [|16|, the energy- 
momentum density tensor is given by: 

< T^,{x) >= 8 J d'p p.F^{x,p) + {f/($) - fv/2 jxf - fp/2 bxb^jg^, 

-^[TrF'ix)^^^ - fvTri^xFixh'Fix)) - fpTr{^xT ■ F{x)^\F{x))]g,, , (5) 
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where F^{x,p) is the isoscalar vector component of the Wigner function and the classical 
value of the a field, obeys the equation ^//g + A + S $^ = ps{x) =<: ilj{x)ip{x) :>. 

The matrix F{x) = J d'^p F{x,p) can be decomposed in components with definite trans- 
formation properties (Clifford algebra), where, for instance, the scalar and the vector 
components of the isoscalar part are related to scalar and current densities as: 

8F{x) = ps{x), 8F'^{x) = j^{x) =<: il){x)'^^il){x) :>. The corresponding components 
of the isovector part are related to isovector densities: h{x) =<: 'ip{x)T^ip{x) :> ( scalar ), 
h^{x) =<: ?/'(x)r37'^'0(x) :> ( vector ). 

The quantities in square brackets of Eq. (5) are the contributions of the exchange terms. 
It is essential to note that Fock terms contain traces of powers of F{x) that naturally bring 
scalar, vector, tensor, pseudoscalar and pseudovector contributions. In particular for the 
case of asymmetric nuclear matter we obtain scalar and vector isovector contributions to 
the EOS, generally associated respectively with 6 and p mesons. From Eq.(6) we obtain 
the energy density for asymmetric nuclear matter that in analogy to the Hartree case can 
be rewritten in the following form: 

e = < Too > = efc./ + e,,^'^ + [/($) + ^fsps' + \~fvPB' + ^/^^^ + ^/(.&o' (6) 

where is the baryon density and = psp — Psn is the corresponding isovector density. 
The 



(27r)3 

are kinetic contributions and 



I £p^|f^^M^ = i [?>E*pB^ + M*pb^-] t = n,p (7) 



1 1 rf$ £^ 3 



1 1 d^^ 1 



5 1 d^ d^^ 3 



are density dependent effective coupling constants. E* = y|pP + M*^ and M* are the 
effective masses, see Eq.s.(9,10). 

Here we explicitly obtain a density dependence arising also in the vector, isovector and 
isoscalar couplings, like in the phenomenological approach of Ref. jl^. We have verified 



that the used approximation leads to a thermodynamically consistent theory 

We remark that, as in non-linear mean- field models, we have in total five parameters. 
As usual, the ones related to isoscalar mesons, fv, fs, A, B are fixed in order to reproduce 
equilibrium properties of symmetric nuclear matter: saturation density po = 0.16 fm~^, 
binding energy E/A = —16 MeV, compressibility modulus Kq = 245 MeV and nucleon 
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effective (or Dirac) mass at po, Mq* = 0.73 M (see Eq.(lO)). This value for Mq* is in the 
range expected from the analysis of elliptic-flow data [18||. The coupling constant /„ can then 



be adjusted in order to get a good value for the symmetry energy at saturation density, but 
now taking into account the contribution to the isovector channel coming from the isoscalar 
mesons through the Fock terms. In our calculations we have a symmetry coefficient of the 
Weiszaecker mass formula = 31.5 MeV. In the Table we list the values obtained for the 
five parameters. 



fv 


fs 


A 

(/m-i) 


B 


fp 


3.998 


9.731 


0.088 


-0.015 


0.6 



Table NLHF parameters from the fit to saturation properties of nuclear 
matter (see the text). 

According to these values, we see that the dominant contributions to the density depen- 
dent coupling functions fs, fv, f'si fvi Eq.(8), come essentialy from the isoscalar a and uj 
mesons. 

We discuss now some results for the EOS of asymmetric nuclear matter. We show the 
comparison between our Non Linear Hartree-Fock (NLHF) present calculations and those 
of the Non Linear Hartree {NLH) model of Ref. [0,0], including the isovector p and 5 
mesons [0], with parameters chosen in order to give the same saturation properties. 

The comparison for potential symmetry energies per nucleon is presented in Fig.l. As 
already stressed, in the NLHF results a large contribution to the symmetry term comes 
from the Fock contributions associated with the a and u mesons, with the corresponding 
four parameters of the theory fitted on properties of symmetric nuclear matter. Therefore 
the inclusion of the Fock terms (solid line, NLHF) can give the correct value of even with 
a relatively small coupling constant for the p meson: fp = 0.6 /m^, close to the free space 
value. We remind that this large exchange contribution to the symmetry energy occurs also 
in QHD without non-linear terms 0. 

We show also NLH calculations including both p and the isovector scalar S mesons 
(dashed line). In this case the coupling constants fp and fg have been chosen in order to 
reproduce at saturation density the same symmetry energy and neutron-proton effective 
mass splitting that we get within our model. We remark that now we need a /p = 2.3 /m^, 
about four times the free space value, and also a relatively strong 6 coupling, fs = 1.4 /m^, 
but still in the range of free space values |19,20]. For reference we show also the result 



of a NLH calculations including only the p contribution (long-dashed line). In order to 
have the same 04 value at saturation density we need a /p = 1.2 frn^, still almost two 
times the free space value. We stress that the inclusion of the 6 contribution in the Hartree 
scheme, necessary for the neutron-proton mass splitting, gives also an attractive term in the 
symmetry energy, see Ref. and so a much stronger p coupling is required in order to 
reproduce the correct 04 coefficient around saturation density. 

In all these relativistic models a quite repulsive density dependence of the symmetry 
term of the EOS is obtained. 

We notice that the density dependence of the symmetry energy that one obtains in the 
complete NLH + p + 6 model is quite different with respect to our results. This is due to 
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the fact that in the NLHF model the couphng functions in the isovector channels (/_^, fy) 
become density dependent. This represents a qualitative new effect of the exchange terms 
in a non-linear QHD model. 

Such density dependence is shown in Figure 2. In particular we stress the behaviour at 
sub-nuclear densities due to the opposite sign of the d^/dps contribution, see Eq(9). This 
implies a "softer" behaviour of the potential symmetry term below saturation density in the 
NHLF case [see the insert in Fig.l]. 

In Fig. 3 we report the density dependence of neutron (bottom) and proton (top) effective 
masses for various asymmetries (/ = {N — Z)/A) as predicted by NLHF (solid lines) and 
NLH + p + 6 (dashed lines). We remind that in the usual Hartree approximation this effect 
is associated with the scalar isovector 6 meson. The Fock terms lead to a behaviour: 

M:,,{NLHF) = M*±f^b + i+^p,-^n), (9) 

where M* is the nucleon effective mass in symmetric nuclear matter, 

M* = M - $ - (/™ + 2f,)ps + ^ipl + Pl)^- (10) 

and 

fm ^ _ 1^ _ fp_ 

^ 2 %dps 2' ^ ' 

Since the coefficient fg^ is positive we get an effect very similar to what expected from the 



contribution of the 5 meson []T9|, dashed lines in Fig. 3. The splitting of proton and neutron 
effective masses influences also the density behaviour of the symmetry energy (Fig.l) and is 
responsible for its rapid increase at high density. On this point we would like to make two 
more remarks: 

i) The splitting is quite small around normal density, so it can be neglected in finite 
nuclei. This difference in n and p effective masses can however be relevant for transport 
properties of asymmetric, dense nuclear matter that can be formed in intermediate energy 
reactions with radioactive beams, naturally apart neutron star properties; 

ii) The proton effective masses are systematically above the neutron ones. This trend. 



also in agreement with 5 calculations |[T9|J21[1 , is just the opposite of what expected from 
Brueckner-Hartree-Fock calculations with realistic A^A^ potentials Although relativis- 
tic and non-relativistic effective masses cannot be directly compared, see ref.s p2|,p|, it 



is interesting to look at the predictions given by Skyrme-like effective forces. Calculations 
based on Skyrme forces are, to some extent, the non-relativistic counterpart of our approach, 
also because exchange contributions are suitably accounted for. Concerning the sphtting of 
neutron and proton effective masses in asymmetric matter, the most recent parametriza- 



tions, SLy-type p4[, of Skyrme forces give the proton effective mass above the neutron one, 
in agreement with our calculations. Previous parametrizations, instead, yield a splitting in 
the opposite direction, but also show unpleasant behaviours in the spin channel (collapse 
of polarized neutron matter, see discussion in [Q). In the insert of Fig. 3 the ratio of the 
splitting, 5M* = M* — M*, to the bare nucleon mass is displayed for the SLy4 force and 
for our approach. 
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The sign of the sphtting depends on the chosen effective interaction. This puzzle can be 
disentangled by a detailed analysis of the transport properties of dense asymmetric nuclear 
matter. 

In conclusion we have shown the evaluation, in a non-perturbative scheme, of Fock term 
contributions in a non linear effective field theory for asymmetric nuclear matter. Very 
reasonable and stimulating results for isospin effects on the nuclear EOS are obtained just 
from such minimal explicit many-body effects. 

We stress again that at variance with non relativistic effective forces all the RMF models 
give a stiff potential symmetry term, i.e. more repulsive with increasing baryon density. 
Moreover the density dependence of the isovector couplings due to the Fock contributions 
leads to a new interesting effect (see the NLHF curve in the insert of Fig.l), a softening of 
the behaviour at sub-nuclear densities. We expect to see related dynamical effects in heavy 
ion collisions at intermediate energies in fragmentation events and collective flows p6| , ^ 



Of course a similar analysis can be performed for spin effects. Moreover a transport 
equation can be consistently derived to be used for the study of dynamical evolution of 
nuclear matter far from normal conditions. 
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FIGURE CAPTIONS 



Figure 1 - Potential symmetry energy per nucleon vs. baryon density (in units of saturation 
density). Solid line: NLHF results. Dashed line: Hartree results with p and 5 mesons 
(NLH). Long-dashed: Hartree results with only p meson. 



Figure 2 - Density dependence of f'g and fy. Each curve is normalized to the value at 
saturation density. 



Figure 3 - Proton (top curves) and neutron (bottom curves) effective masses vs. baryon 
density for various charge asymmetries. / = = : long dashed line. / = 0.8 : solid 



lines NLHF] dashed lines NLH. In the insert, the relative splitting of neutron and proton 
effective masses. Solid line: NLHF results, circles: SLy4 |^| non-relativistic results. 



9 




o 



o 

bo 



to 



o 



CO 





1 1 1 1 1 1 1 1 1 1 1 1 1 




1 


1 




-A- O 








- 


+ 






- 


_ 


-A- O 






_ 


— 


+ o 






— 




-A-n 
O 


















U ^ 








— 


o + 
o + 






— 














— 


n -6- 






— 





o + 
o + 









- 


+ 






- 




o + 
o + 










o + 










o 

o + 


o 








o + 


<1 








o + 




o 






o + 


o 


1 — ' 






o + 


o 


p 






0+ 










o 








1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 


1 


1 



